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ABSTRACT 

We report high precision spectrophotometric observations of four planetary transits 
of HD 209458, in the region of the sodium resonance doublet at 589.3 nm. We find that 
the photometric dimming during transit in a bandpass centered on the sodium feature 
is deeper by (2.32 it 0.57) x 10~^ relative to simultaneous observations of the transit 
in adjacent bands. We interpret this additional dimming as absorption from sodium 
in the planetary atmosphere, as recently predicted from several theoretical modeling 
efforts. Our model for a cloudless planetary atmosphere with a solar abundance of 
sodium in atomic form predicts more sodium absorption than we observe. There are 
several possibilities that may account for this reduced amplitude, including reaction of 
atomic sodium into molecular gases and/or condensates, photoionization of sodium by 
the stellar flux, a low primordial abundance of sodium, or the presence of clouds high 
in the atmosphere. 
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1. Introduction 

Since the discovery of planetary transits in the hght curve of HD 209458 (Charbonneau et 
al. 2000; Henry et al. 2000; Mazeh et al. 2000), this star has been the subject of intensive study: 
Multicolor observations (Jha et al. 2000; Deeg, Garrido, k, Claret 2001) have confirmed the expected 
color-dependence of the light curve Awe to the stellar limb-darkening. Radial-velocity monitoring 
during transit (Queloz et al. 2000; Bundy & Marcy 2000) has yielded variations in excess of the 
orbital motion, due to the occultation of the rotating stellar limb by the planet. Very-high-prccision 
photometry (Brown et al. 2001, hereafter BOl) permitted an improved estimate of the planetary 
and stellar radii, orbital inclination, and stellar limb-darkening, as well as a search for planetary 
satellites and circumplanetary rings. 

Transiting extrasolar planets present a unique opportunity for us to learn about the atmo- 
spheres of the objects. Wavelength-dependent variations in the height at which the planet becomes 
opaque to tangential rays will result in wavelength-dependent changes in the ratio of spectra taken 
in and out of transit. Several groups (Seager & Sassclov 2000; Brown 2001; Hubbard et al. 2001) 
have pursued theoretical explorations of this effect for a variety of model planetary atmospheres. 
These studies have demonstrated that clouds, varying temperature structure and chemical com- 
position, and even atmospheric winds, all produce variations that would be observed, should the 
requisite precision be achieved. Based on these calculations, the expected variations could be as 
large as 0.1% relative to the stellar continuum. 

A previous search by Bundy & Marcy (2000) for variations between spectra of HD 209458 
observed in and out of transit (using extant spectra gathered for radial velocity measurements) was 
sensitive to features with a relative intensity greater than 1-2% of the stellar continuum. Moutou 
et al. (2001) conducted a similar study, and achieved a detection threshold of 1%. Although the 
precision achieved by both these studies was insufficient to address any reasonable model of the 
planetary atmosphere, Moutou et al. (2001) placed limits on the models of the planetary exosphere. 
Searches for absorption due to a planetary exosphere have also been conducted for 51 Peg during 
the times of inferior conjunction (Coustenis et al. 1998; Rauer et al. 2000). 

2. Observations and Data Analysis 

We obtained 684 spectra^ of HD 209458 with the HST STIS spectrograph, spanning the times 
of four planetary transits, on UT 2000 April 25, April 28-29, May 5-6, and May 12-13. The 
primary science goal of this project was to improve the estimate of the planetary and stellar radii, 
orbital inclination, and stellar limb-darkening, as well as conduct a search for planetary satellites 
and circumplanetary rings; we presented these results in BOl. In selecting the wavelength range 
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for these observations, the desire to maximize the number of detected photons led us to consider 
regions near 600 nm, where the combination of the instrumental sensitivity and the stellar flux 
would be optimized. 

The secondary science goal was to pursue the prediction by Seager & Sasselov (2000) (and, 
later, Brown 2001; Hubbard ct al. 2001) of a strong spectroscopic feature at 589.3 nm due to 
absorption from sodium in the planetary atmosphere. Thus we chose to observe the wavelength 
region 581.3 < A < 638.2 nm, with a medium resolution of i? = A/AA = 5540, corresponding to a 
resolution element of 0.11 nm. 

The details of the data acquisition and analysis are presented in BOl, and we refer the reader to 
that publication. In summary, the data reduction consisted of (1) recalibrating the two-dimensional 
CCD images, (2) removing cosmic-ray events, (3) extracting one-dimensional spectra, (4) summing 
the detected counts over wavelength to yield a photometric index, and (5) correcting the resulting 
photometric time series for variations that depend on the phase of the HST orbit and for variations 
between visits. The only difference in these procedures between BOl and the current work is in 
step (4): Previously, we summed the spectra either (a) over the entire available wavelength range, 
or (b) over the blue and red halves of the available range. In the present paper, we restrict greatly 
the wavelength span over which we perform the integration. 

As described in BOl, observations of the first transit (UT 2000 April 25) were partially com- 
promised by a database error in the location of the position of the spectrum on the detector. The 
result was that the spectrum was not entirely contained within the CCD subarray. In the subse- 
quent data analysis, we ignored these data. Furthermore, as described in BOl, the first orbit (of 
5) for each visit shows photometric variability in excess of that achieved for the remaining orbits. 
These variations are presumably due to the spacecraft; we omit these data as well. This leaves 417 
spectra of the 684 acquired. 

Since we do not know the precise width of the feature we seek, we select three bands of varying 
width, each centered on the sodium feature. We refer to these bands as "narrow" (n), "medium" 
(m) , and "wide" {w) . The n band is the smallest wavelength range that still encompasses the stellar 
sodium lines; the w band is the widest wavelength range that permits an adjacent calibrating band 
to the blue. The m, band is roughly intermediate (by ratio) between these two extremes; it is ~ 1/3 
the range of the w band, and ~ 3 times the range of the n band. We further define, for each of 
these, a "blue" (6) and a "red" (r) band, which bracket the "center" (c) band. The names and 
ranges of these nine bands are given in Table 1, and displayed in Figure 1. For each of these 9 
bands, we produce a photometric time series by the procedure described above. The photometric 
index at a time t is then identified by a letter indicating the width, and a subscript indicating the 
position (e.g. ^^Wr{t)" indicates the photometric index in the wide band, red side, at time t). Each 
of these is a normalized time series, with a value of unity when averaged over the out-of-transit 
observations; the minimum values near the transit centers are approximately 0.984. 

We denote the time of the center of the photometric transit by Tc. In what follows, we 
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consider the in-transit observations, (i.e. those that occur between second and third contacts, 
\t — Tc\ < 66.111 min), which we denote by and the out-of-transit observations (those that 
occur before first contact, or after fourth contact, \t — Tc\ > 92.125 min), which we denote by tout- 
There are 171 in-transit observations, and 207 out-of-transit observations. We ignore the small 
fraction (7%) of observations that occur during ingress or egress. 



2.1. Stellar Limb-Darkening 

One potential source of color-dependent variation in the transit shape is stellar limb-darkening 
(see Figure 6 in BOX; also Jha et al. 2000; Deeg, Garrido, & Claret 2001). In order to investigate 
this possibility, we produce (for each width) the difference of the red and blue bands: 

Udit) = nb{t) - nr(t) 

md{t) = mb{t) - mr{t) (1) 

Wd{t) = Wb{t) - Wr{t). 

The observed standard deviations in these time series (measured over the out-of-transit observa- 
tions) are cr[nd{tout)] = 3.3 x 10"^, (T[md{tout)] = 3.6 x 10~^, and a[wd{tout)] = 5.4 x 10"^. These 
values match the predictions based on photon noise. 

In order to check for changes in the transit depth due to stellar limb-darkening, we then 
calculate, for each of the time series above, the difference in the mean of the relative flux, between 
the in-transit observations and the out-of-transit observations: 

And = ridjUn) - nditput) = (-2.0 ± 2.5) x 10"^ 

Amd = mdjtin) - mditout) = (+0.1 ± 2.7) x 10"^ (2) 

Awd = WdiUn) - Wditout) = (-3.1 ± 4.1) X 10-5. 

Overbars indicate averages over time, and each quoted error is the estimated standard deviation 
of the mean. All of these values are consistent with no variation. Thus, we have no evidence for 
color-dependent limb-darkening between these bands, which span about 15 nm. 

We compared these results to predictions based on a theoretical model of the stellar surface 
brightness. The model was produced by R. Kurucz (personal communication, 2000). We denote it 
by S{X, fi), where fi is the cosine of the angle between the line of sight and the normal to the local 
stellar surface. The stellar parameters were taken to be the best-fit values as found by Mazeh et 
al. (2000): T^ff = 6000 K, logg = 4.25, and [Fe/H] = 0.00. The model is evaluated at a spectral 
resolution of R = 2 000 000 (greatly in excess of the STIS resolution, R = 5540), and at 17 values 
of /i G [0, 1] 

The stellar model allows us to produce a theoretical transit curve over a chosen bandpass, 
which includes the effects of stellar limb-darkening. We do not force the limb-darkening to fit 
a parameterized model. Rather, we interpolate between the values of /x on which the model of 
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S{X,ijl) is initially calculated. For each A, we produce the theoretical transit curve, T{X,t), for all 
t at which we have data. The curve is calculated by integrating over the unocculted portion of the 
stellar disk as described in Charbonneau et al. (2000) and Sackctt (1999). Wc assume the best-fit 
values from BOl for the stellar radius (Rs = 1.146 Rq), planetary radius (Rp = 1.347 -Rjup, where 
Rjup = 71492 km is the equatorial radius of Jupiter at a pressure of 1 bar; Cox (2000)), orbital 
period {P = 3.52474 days), orbital inclination {i = 86?6) and semi-major axis (a = 0.0468 AU), 
assuming a value for the stellar mass of Mg = 1.1 Mq (Mazeh et al. 2000). 

We note that the STIS data have an effective continuum that is tilted relative to the theoretical 
values, due to the wavelength dependent sensitivity of the instrument. We wish to emulate the 
data, where each A is weighted by the observed number of photons. To do so, we fit a low-order 
polynomial (in A) to the model disk-integrated spectrum, and a second low-order polynomial to 
the STIS data. We subsequently divide the calculated T(X,t) for each t by the ratio of these 
polynomials. 

Next, we integrate T(X,t) in A over each of the band passes in Table 1. We subsequently 
difference the model light curve as calculated for the blue and red bands (emulating what we did 
with the data in eq. [1]), and then evaluate the mean of the difference (as in eq. [2]): 

Ana (theory) = -0.44 x 10"^ 

Anid (theory) = -0.58 x 10"^ (3) 
Awd (theory) = -0.85 x 10"^. 

These results can be compared directly with equation (2). Each of these is much smaller than 
the observational precision we achieved, and thus consistent with our observational result of no 
significant offset. 

Furthermore, the results stated in equation (2) demonstrate that the data remain photon- 
noise limited at these high levels of precision: We have taken the mean of two large groups of data 
(roughly 200 observations apiece), and found that the difference is consistent with photon- noise 
limited photometry, with a typical precision 3 x 10~^. 



2.2. The Sodium Band 

In order to search for variations in the sodium band relative to the adjacent bands, we produce 
(for each width) the mean light curve of the blue and red bands, and difference this from the light 
curve for the center band: 

nNa{t) = ndt) - [nb{t) + nr{t)] / 2 

rriNait) = mc{t) - [mb{t) + mr{t)] / 2 (4) 

WNait) = Wc{t) - [Wb{t) + Wr{t)] / 2. 

This linear combination removes the variations due to the color dependence of the limb-darkening of 
the stellar continuum. As we found above (§ 2.1), this effect is very small (eq. [3]). We consider the 
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effect of the deviations in the hmb-darkening that occur in the cores of stellar absorption features in 
§4.1, and show that this effect is also negligible. The three time scries in equation (4) are plotted as 
a function of absolute value of the time from the center of transit in Figure 2. The observed standard 
deviations in these time series over the out-of-transit observations are ulnNaitout)] = 5.5 x 10~^, 
c"["T'iVa(*o«t)] = 3.7 X 10~^, and (T[w]s[a{tout)] = 3.5 x 10"^. These achieved values match the 
predictions of the photon-noise-limited precision. 

In order to look for changes in the transit depth in the sodium-band relative to the adjacent 
bands, we calculate, for each of the time series above, the difference in the mean as observed in and 
out of transit: 

Aniva = riNaitin) - nNajtout) = (-23.2 ± 5.7) x 10-5 

AmATa = rnNajtin) - mNajtout) = (-13.1 ± 3.8) X 10"^ (5) 
AwNa = WNaiUn) - WNa{tout) = (- 3.1 ± 3.6) X 10"^. 

As before, the errors shown are the 1-a errors of the mean. The results indicate that we have 
detected a deeper transit in the sodium band for the narrow and medium bandwidths, with a 
significance of 4.1 a and 3.4 a, respectively. We find no significant offset for the wide bandwidth. 

Furthermore, wc calculate the three quantities in equation (5) for each of the three transits 
separately. These results confirm the conclusion that there is a deeper signal in the narrow and 
medium bands. The variations between visits for the same band are not significant given the 
precision of the data. 

In Figure 3, we plot histograms of the in-transit data and the out-of-transit data, for the narrow 
band. This plot shows that each of these sets of points appears to be drawn from a normal distribu- 
tion with a standard deviation as predicted by photon-counting arguments. However, the in-transit 
points scatter about a mean that is significantly offset from that defined by the out-of-transit points. 
This indicates that we have not skewed these distributions by our analysis procedures: Rather, the 
in-transit observations appear to represent a shifted version of the out-of-transit observations. 

We bin njva(i) in time, and plot these results in Figure 4. These further illustrate that we 
have observed a deeper transit in the sodium band. 

We investigate the possibility that the observed decrement is due to non-linearity in the STIS 
CCD: Since the observed value of Ara^va is —2.32 x 10~^ for a change in mean intensity inside 
to outside transit of —1.6 x 10"^, a non-linearity of ~ 1% across this range would be required. 
Gilliland, Goudfrooij, &; Kimble (1999) derive an upper limit that is an order of magnitude lower 
than this. We further test this effect directly by selecting 18 strong stellar absorption features in 
our observed wavelength range. For each spectral feature, we define a set of band passes with the 
same wavelength range as those listed in Table 1, but now centered on the spectral line. We derive 
An and Am for each of these as in equation (5). The results are shown in Figure 5. We find no 
correlation of either An or Am with spectral line depth. These tests also confirm our evaluation 
of the precision we achieve in equation (5): We find that the values of An or Am are normally 
distributed, with a mean consistent with 0, and a standard deviation as we found above (eq. [5]). 
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3. Comparison with Theoretical Predictions 

As described above, the quantity that we observe is the difference in the transit depth in a 
band centered on the Na D hnes to the average of two flanking bands, as a function of time from the 
center of transit: We denoted these by njva(i)) i^Nait), and WNait): above. We wish to compare 
these results to models of the planetary atmosphere. Several steps arc required to transform the 
model predictions into the same quantity as the observable, and we describe these steps here. 

We first produce several model calculations of the change in the effective planetary radius as 
a function of wavelength, as prescribed by Brown (2001). For all the models, we use the best-fit 
values from BOX for Rp, Rg, i, a, Mg as given in § 2.1. We further specify the planetary mass 
{Mp = 0.69 Mjup) and the stellar effective temperature {Tg = 6000 K), both as given by Mazeh et 
al. (2000). We set the equatorial rotational velocity of the planet to Vrot = 2.0 kms"-*^. This is the 
value implied by the measured planetary radius, under the assumption that the planet is tidally 
locked (i.e. the rotational period is the same as the orbital period). Each model includes the effect 
of photoionization as described in Brown (2001), which reduces the core strength of the Na D 
lines, but does not significantly change the wings. 

The fiducial model has a cloud deck with cloud tops at a pressure of .0368 ( = 0.1 / e) bar, 
and solar metallicity (the measured value for HD 209458; Mazeh et al. 2000). We also consider a 
number of variants to the fiducial model, representing changes to either the metal abundance, or 
the cloud height. These variants are listed in Table 2. For each model (denoted by the index m), 
we have the theoretical planetary radius Rp{\,m), where A denotes the wavelength. For each of 
these, we renormalize Rp{X,m) so that the average value of the planetary radius over the entire 
STIS bandpass is the best-fit value determined by BOl (1.347 -Rjup)- These renormalized functions 
Rp{X,m) are shown in Figure 6. 

We use the theoretical model of the stellar surface brightness [S{X,n); described in § 2.1] to 
include the effects of the wavelength-dependent limb-darkening of the star. This approach allows us 
to account for both the continuum limb-darkening, and the deviations that occur from this in the 
cores of absorption lines; these can be significant in the case of the sodium features. In general, the 
limb-darkening is less pronounced in the cores of the sodium lines than in the adjacent continuum, 
and thus for a planetary radius that is independent of wavelength, the transit light curve near 
Tc is deeper and more rounded as observed in the continuum than as observed in the core of an 
absorption line. 

We proceed as described in § 2.1, with the exception that Rp{\,m) is now a function of 
wavelength: For each A, wc produce the theoretical transit curve, T{\,t,m). We bin these results 
in the wavelength bands as described in § 2, and difference the values in the band centered on the 
sodium lines from the average of the flanking bands, as prescribed by equation (4). The result is 
a theoretical time series that closely resembles the observed STIS time series, with the effect of a 
wavelength-dependent planetary radius included. These model time series are shown in Figure 7, 
and may be compared directly with the data shown in Figures 2 &; 4. 
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We find that most of the models hsted in Table 1 produce differential transits that are signif- 
icantly deeper than our observational result. Only models n4 or c3, which represent, respectively, 
extreme values of high cloud height or depleted atomic sodium abundance, produce transits of 
approximately the correct depth. 



We have detected a significantly deeper transit in the sodium band relative to the adjacent 
bands. We interpret the signal as due to absorption by sodium in the planetary atmosphere. We 
are encouraged in this interpretation by the following two considerations: First, a sodium feature 
in the transmission spectrum of HD 209458 was predicted unanimously by the current modeling 
efforts (Seager & Sasselov 2000; Brown 2001; Hubbard et al. 2001). Second, observed spectra of L- 
type objects (e.g. Kirkpatrick et al. 1999), T-dwarfs (e.g. Burgasser et al. 1999) and, in particular, 
the brown dwarf Gl 229 B (e.g. Oppenheimer et al. 1998), show very strong absorption from 
alkali metal lines (Burrows, Marley, &; Sharp 2000). These objects span a temperature range of 
900-2000 K, which includes the equilibrium temperature of HD 209458 b, Teff = 1430 (1 - A)i K 
(where A indicates the Bond albedo). (It should be noted that these objects are at roughly the same 
temperature as HD 209458 b for very different reasons: The planet is heated by its proximity to the 
parent star, whereas the temperatures for the other objects are set by their ongoing contraction.) 

In the following discussion, we first consider alternate interpretations of our result, then discuss 
the constraints that we can place on the planetary atmosphere, and finish by looking ahead to near- 
future complementary observations. 



We must take care that the observed decrement is not the result of the distinctive limb- 
darkening exhibited by the sodium line relative to that of the adjacent continuum. In order to 
quantify the amplitude of this effect, we generated model transit curves (for a planet of constant 
radius) as described in § 2.1, and integrated these in A over the nine band passes listed in Table 1. 
We then recreated the linear combination of these results as given in eq (4). We found that the 
predicted offsets were: 



4. 



Discussion 



4.1. Alternate Explanations 



AnjVa (theory; constant i?p) 
^iTT-Na (theory; constant Rp) 
AtUATa (theory; constant Rp) 



+1.52 X 10-5 
+0.39 X 10-5 
+0A7 X 10-5. 



(6) 



These values are much smaller than the signal we detect. Moreover, they are of the opposite sign: 
Since the net decrement toward the limb is less in the sodium line than in the adjacent continuum, 
the transit as observed in a band centered on the sodium lines would be expected to be less deep near 
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Tc than that observed in a band in the adjacent spectral region. We note that the above quantities 
should be subtracted from the observed values of An^ra, Arriba, and Awj^a given in equation (5) for 
comparison with future modeling efforts that do not explicitly include the wavclcngth-dcpcndcnt 
limb-darkening. For the purposes of the present paper, our models explicitly account for this effect 
by making use of the model stellar surface brightness, and thus no adjustment is needed. 

A second concern is that the star might appear smaller when observed at the wavelength of 
the sodium lines, producing a deeper transit for the same-sized planet. In fact, we expect the 
opposite to be true: The solar limb should be slightly larger in the core of the Na D lines because 
of the great opacity in these lines. Model calculations (E. H. Avrett, personal communication) 
yield an intensity-weighted solar radius increase, after averaging over even the narrowest bandpass 
containing the Na D lines, of less than 10 km. This will have a completely negligible effect on the 
observed transit depth in the bandpass containing the Na D lines. 

4.2. Constraints on the Planetary Atmosphere 

Our fiducial model (si), which has cloud tops at 0.04 bar and a solar abundance of sodium in 
atomic form, predicts values for An Na and Am Na that are ~3 times deeper than we observe. This 
conclusion awaits confirmation from other more detailed models, such as those performed by Seager 
& Sasselov (2000) and Hubbard et al. (2001). Here we consider several possible physical effects 
that may contribute to the reduced sodium absorption relative to the predictions of our model: 

It may be that a significant fraction of atomic sodium has combined into molecules, which may 
either be present as gaseous species, or sequestered from the atmosphere as condensates. Chemical 
equilibrium models for substellar and planetary atmospheres (Fegley & Loddcrs 1994, 1996; Burrows 
&: Sharp 1999; Lodders 1999) indicate that Na2S, NaCl, NaOH, NaH, and more complicated species 
(such as NaAlSiaOs) may be present in the planetary atmosphere. In particular, Lodders (1999) 
finds that, at these temperatures, the competition for sodium is predominantly between Na2S 
(solid) and monatomic Na (gas). At temperatures above the Na2S-condensation line, sodium is 
predominantly in atomic form, and gases such as NaCl, NaOH, and NaH never become the most 
abundant sodium-bearing species (K. Lodders, personal communcation) . If the disparity between 
the observed values for AriNa and AniNa and the predictions for these quantities based on our 
fiducial model (si), is due entirely to this effect (and not, for example, mitigated by the effect of 
high cloud decks), then this would indicate that less than 1 % of the sodium is present in atomic 
form. This statement is based on the fact that even our model c3 fails to reproduce the observed 
values. 

Our current models (Brown 2001) account for the ionization of sodium by the large stellar 
flux incident upon the planetary atmosphere. However, it is possible that the effect is larger than 
predicted by our calculations: This would require a sink for free electrons, the effect of which would 
be to lower the recombination rate. 
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Perhaps HD 209458 began with a depleted sodium abundance, or, more generally, a depleted 
metal abundance. The current picture of gas-giant planet formation predicts that planets should 
present a metallicity at or above the value observed in the parent star; Jupiter and Saturn are 
certainly metal-rich relative to the Sun (see discussion in Fegley & Lodders 1994). Metallicity 
appears to play an important role in the formation of planets: Parent stars of close-in planets 
appear to be metal rich (Gonzalez 1997, 2000) relative to their counterparts in the local solar 
neighborhood (Queloz et al. 2000b; Butler et al. 2000). Furthermore, the observed lack of close-in 
giant planets in the globular cluster 47 Tuc (Gilliland et al. 2000) may be due, in part, to a reduced 
planet-formation rate due to the low-metallicity of the globular cluster. The relationship linking 
metallicity and planet formation is not clear. 

Clouds provide a natural means of reducing the sodium absorption, by creating a hard edge to 
the atmosphere, and thereby reducing the effective area of the atmosphere viewed in transmission. 
If the difference between the observed values of An^a and Arriba and the predictions for these 
quantities based on our fiducial model (si) is due entirely to clouds, then this would require a very 
high cloud deck, with cloud tops above 0.4 mbar (model n4). Ackerman &; Marley (2001) have 
recently produced algorithms to model cloud formation in the atmospheres of extrasolar planets. 
The effect of the stellar UV flux on the atmospheric chemistry has not been explicitly included. 
We urge investigations of the production of condensates, such as photochemical hazes, by this 
mechanism. Such photochemical hazes are present in the atmospheres of the gas giants of the solar 
system, and produce large effects upon the observed reflectance spectra of these bodies (Karkoschka 
1994) . The stellar flux incident upon the atmosphere of HD 209458 b is roughly 20,000 times greater 
than that upon Jupiter. 



4.3. Future Observations 

In the present work, we have demonstrated that by using existing instruments, it is possible 
to begin the investigation of the atmospheres of planets orbiting other stars. Wide-field surveys 
conducted by dedicated, small telescopes (e.g. Borucki et al. 2001; Brown & Charbonneau 2000)^ 
promise to detect numerous close-in, transiting planets over the next few years. The parent stars 
of these planets should be sufficiently bright to permit similar studies of the planetary atmosphere. 

Based on the results presented here, we believe the observed value of sodium absorption from 
the planetary atmosphere may have resulted from a high cloud deck, a low atomic sodium abun- 
dance, or a combination of both. Observations of the transmission spectrum over a broader wave- 
length range should allow us to distinguish between these two broad categories of models (Figure 8). 
If the effect is predominantly due to clouds, the transmission spectrum from 300-1000 nm will show 
little variation with wavelength, other than additional alkali metal lines. If, on the other hand, the 



'^See also http://www.hao.ucar.edu/public/research/staxe/stare.html, maintained by T. Brown and D. Kolinski. 
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clouds are very deep or not present, then the transmission spectrum may show deep features red- 
ward of 500 nm due to water, and absorption due to Rayleigh scattering blueward of 500 nm. These 
arc only two extremes of the many possibilities; the actual case may share characteristics of both 
of these effects. Further diagnostics could be provided by near-IR spectroscopy, which could detect 
features due to molecular species (such as H2O, CO, and CH4) that are expected to be present 
in planetary transmission spectra (Seager &; Sasselov 2000; Brown 2001; Hubbard et al. 2001). 
Suitable observations in the 1-5 /xm region would allow inferences about the atmosphere's ther- 
mal structure and composition, permitting a more definite interpretation of the sodium signature. 
Brown et al. (2000) have presented preliminary near-IR spectra of HD 209458 in and out of transit. 
These spectra demonstrate that it should be possible to search for near-IR atmospheric signatures 
with existing instruments, such as the NIRSPEC spectrograph (McLean et al. 1998) on the Keck II 
telescope. 

Observations of the reflected light spectrum (wavelength-dependent geometric albedo) as pur- 
sued for the r Boo system (Charbonneau et al. 1999; Collier Cameron et al. 1999) would complement 
the results of transmission spectroscopy. Predictions for the albedo as a function of wavelength vary 
by several orders of magnitude (Seager & Sasselov 1998; Marley et al. 1999; Goukcnlcuque et al. 
2000; Sudarsky, Burrows, & Pinto 2000). These observations should be facilitated in the case of a 
transiting system (Charbonneau & Noyes 2000). If possible, the observation of the phase function 
(variation in the reflected-light with orbital phase) would constrain greatly the nature of particu- 
lates in the planetary atmosphere (Seager, Whitney, & Sasselov 2000). Observations of molecular 
absorption features from thermal emission spectra, such as sought in the recent study of the r Boo 
system (Wiedemann, Deming, & Bjorakcr 2001), would clearly complement these studies. Finally, 
infrared photometry of the secondary eclipse would provide a direct measurement of the effective 
temperature of the planetary atmosphere, should the requisite precision be achieved. 
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Table 1. Wavelength Bands 



width 


blue (nm) 


center (nm) 




red (nm) 


narrow 


rib-. 581.8-588.7 


ric: 588.7-589.9 


rir. 


589.9-596.8 


medium 


mbi 581.8-587.4 


mci 587.4-591.2 


rrir 


: 591.2-596.8 


wide 


iri,: 581.8 584.3 


11',: 584.3 594.3 


w,.: 


594.3 596.8 



Tabic 2. Model Atmosphere Parameters 



identifier 


cloud tops 
(bar) 


metal abundance 
(relative to solar) 


si 


0.0368 


1.00 


nl 


3.68 


1.00 


n2 


0.368 
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n3 


0.00368 


1.00 


n4 
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1.00 
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0.0368 
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c2 
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0.0368 
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Fig. 1. — Shown is a portion of a STIS spectrum of HD 209458, centered on the Na D hnes. The 

vertical axis is the number of detected photoelcctrons per wavelength pixel after integrating over 
17 pixels in cross-dispersion. The vertical lines indicate the band passes over which we integrate 
the spectrum to produce the photometric time series. The band ric is the set of pixels between the 
two dotted lines; nic is the set between the dashed lines; Wc is the set between the dot-dashed lines. 
The corresponding blue bands {rib, frib, and Wb) are the set of pixels between the left solid line and 
the left boundary of the center band. Similarly, the corresponding red bands (n^, rur, and Wr) are 
the set of pixels between the right edge of the center band and the right solid line. 
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Fig. 2. — (Previous page) The unbinned time series njva, mNa, and WNa are plotted in the top, 
center, and bottom panels, respectively, as a function of absolute time from the center of transit. 
The mean of the in-transit values of njva and iriNa are both significantly offset below 0. 
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Fig. 3. — The dark solid line is a histogram of the out-of-transit values of njva (the data shown in 
the upper panel of Figure 2). The dashed dark curve is a Gaussian distribution with a mean of 0, 
o = 5.5 X 10^^ as prescribed from photon-noise-prcdictions, and normalized to the same area. The 
gray solid line is a histogram of the in-transit values of riMa^ renormalized to the same number of 
observations as the dark solid line. The dashed gray curve is a Gaussian distribution with the same 
(7 as the dark curve, but with a mean of —2.32 x 10"^ (the observed value). The gray histogram is 
significantly offset from 0. 
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Fig. 4. — The upper panel shows the unbinned time series n^a (the data shown in the upper panel 
of Figure 2). The lower panel shows these data binned in time (each point is the median value in 
each bin). There are 10 bins, with roughly equal numbers of observations per bin (42). The error 
bars indicate the estimated standard deviation of the median. The solid curve is a model for the 
difference of two transit curves (described in § 3) , scaled to the observed offset in the mean during 
transit, /^nna = -2.32 x 10"^. 
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Fig. 5. — (Previous page) The upper panel shows a typical spectrum of HD 209458, over the full 
spectroscopic range that we observed. We selected 18 strong features in addition to the Na D lines; 
these arc indicated by vertical dashed lines. The middle panel shows the results for the narrow 
(n) band pass. The black diamond is the value of AnA^a- The wavelength ranges of the bands rib, 
He, and rir are shown as the horizontal gray bars. We then repeated the analysis, i.e. same band 
widths, but now centered over each of the 18 features that are indicated. The resulting values are 
shown as white diamonds. The lower panel shows the corresponding results for the middle (m) 
band pass. In each case, the white diamond points appear to be normally-distributed, with a mean 
of 0, and a standard deviation as prescribed by photon statistics. In particular, the residuals do not 
correlate with the depth of the spectral feature. Some correlations between adjacent white points 
are evident, as expected, since the band passes are wider than the separation between points. In 
the middle and lower panels, the black diamond point is inconsistent with no variation at the level 
of 4.1 a and 3.4 a, respectively. Note that the y-axis in the lower panel has been scaled so that the 
error bars appear the same size as in the middle panel. 
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Fig. 6. — The effective planetary radius is shown as a function of wavelength for a variety of 
atmospheric models; the vertical lines indicate the bandwidths, and are the same lines as those 
in Figure 1. In both panels, the dark line shows the planetary radius for the fiducial model (si). 
In the upper panel, we show the variation with cloud height; the upper and lower gray curves are 
models n3 and n4, respectively. (The results for models nl and n2 are indistinguishable on this 
plot from that for si). In the lower panel, we show the variation with sodium abundance; the gray 
curves are, from top to bottom, models cl, c2, and c3. 
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Fig. 7. — (Previous page) The solid curve in the top panel shows the model curve for riNa for the 
fiducial model (si). The solid curves in the middle and lower panels are the results for mj\fa and 
wjsfa, respectively, for this same model. The dashed lines show the variations due to cloud height: 
Within each panel, the upper dashed curve is for model n4, and the lower dashed curve is for model 
n3. (The results for models nl and n2 are indistinguishable on this plot from that for si). The 
solid gray lines show the variations due to atomic sodium abundance: Within each panel, the gray 
curves correspond (from top to bottom) to models c3, c2, and cl. These curves may be compared 
directly with the data presented in Figures 2 and 4. 
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Fig. 8. — Shown are three model transmission spectra for HD 209458 over a wavelength range that 

exceeds that of our current data set. The fiducial model (si) with a solar abundance of sodium, and 
cloud tops at 0.0368 bar, is shown as a dotted line. This model is excluded by the data, which do 
not permit such a deep sodium feature. One possibility is that very high clouds reduce the depth 
of the sodium absorption feature; such a model is shown as the solid line. Another possibility is 
that the abundance of atomic sodium in the planetary atmosphere is reduced greatly from the solar 
value; this is shown as the gray line. Observations of the transmission spectrum over the range 
shown here should allow for us to distinguish between these two broad categories of models. 



